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CALCULATING  THE  MAGNETIC  FIELD  FROM  A  SADDLE-SHAPED  COIL  WITHOUT 
A  FERROMAGNETIC  CIRCUIT 


L.  P.  Gnedin,  Dr.  Tech.  Sci . ;  G.  G.  Borzov,  Engr . 


Mathematical  formulae  for  the  components  of  induction 
of  the  magnetic  field  and  the  magnetic  flux  from  a  coil 
with  saddle-shaped  side  sections  in  a  system  without  a 
magnetic  circuit  are  given. 


The  use  of  superconductors  to  generate  high-strength  magnetic 
fields  in  the  excitation  systems  of  MHD  converters  and  electrical 
machines  makes  it  possible  to  eliminate  the  magnetic  circuit  in  these 
devices  [ 1 ] . 

The  most  popular  configuration  of  the  excitation  coil  of  an  elec¬ 
trical  machine  can  be  composed  of  two  elements:  a  rectilinear  side 
part,  and  an  arc-shaped  front  section  with  radius  R  (Fig.  1)  -  the 
so-called  C  type  [2].  The  problem  of  calculating  the  induction  of  the 
magnetic  field  from  the  arc-shaped  front  section  was  solved  using  the 
Biot-Savart  law  in  [3]*  in  which  the  approximate  formulae  proposed  in 
[3]  are  only  satisfied  for  small  distances  from  the  Z-axis. 

The  effect  of  the  front  sections  was  not  considered  in  [4,  5]. 

In  short  coils,  e.g.,  in  the  exci tat  ion  coils  of  a  DC  collector  or  in 
MHD  converters  with  a  short  channel,  the  component  of  induction  from 
the  front  section  can  exceed  the  induction  from  the  side  part. 
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Fig.  1.  Calculating  the 
magnetic  field  of  a  sad¬ 
dle-shaped  coil. 


In  order  to  calculate  the  magnetic  field 
of  an  electrical  machine,  it  is  necessary  to 
know  both  the  induction  at  any  point  in  space, 
e.g.,  for  determining  the  losses  in  the  arma¬ 
ture  circuits,  and  in  the  main  working  flux 
through  the  polar  surface  of  the  armature. 

During  an  engineering  calculation  of  an 
electrical  machine  without  usin^  an  electron¬ 
ic  computer  or  for  checking  a  program  which 
has  been  written,  it  is  desirable  to  express 
the  induction  and  flow  using  elementary 
functions.  This  problem  is  considered  in 
the  proposed  article. 


We  will  define  the  magnetic  induction  as  the  vector  potential  A 
in  cylindrical  coordinates: 

;♦  dz,  I 

o  _  '  1  I  f 

*  ?♦  d?,  ?.  ^ 

dA,  dA, 

where  are  the  coordinates  of  the  current  point. 

•The  origin  of  coordinates  coincides  with  the  rotational  axis  of 
the  armature,  and  the  counting  axis  of  the  angles  passes  through  the 
center  of  the  polar  surface  of  the  armature. 

We  will  make  the  following  assumptions. 

1.  We  will  replace  the  cross  section  of  the  coil  with  an  infinite¬ 
ly  thin  conductor  at  the  geometric  center  of  the  cross  section  with 
current  i  equivalent  to  the  current  of  the  entire  cross  section  of 

the  coil. 

2.  We  will  consider  the  plane  in  which  the  front  part  is  located 
to  be  perpendicular  to  the  side  parts  of  the  coil. 


2 


INDUCTION  COMPONENTS  FROM  FRONT  PART  OF  EXCITATION  COIL 
(from  circuit  2  in  Fig.  1) 


a)  The  radial  induction  component  is  The  component  Az  =  0 

in  expression  (1)  for  B.,  (assumption  2).  For  the  direction  of  the 
currents  used  in  Fig.  1,  the  component 


cos  (y,  —  f)  df 

V  fo  +  -  2Pf,  cos  (f,  -  ?)  4-  (i/2  +  z,)«  ’ 


(2) 


where  <p«  is  the  angle  of  placement  of  the  side  part  of  the  excitation 
coil;  T  is  the  angular  position  of  the  current  element;  L  is  the 
length  of  the  excitation  coil. 


Substituting 
we  obtain: 

where 


(2)  in  expression  and  setting 


.a.i*Mi/2-f 'Z.) 
iinsl/T 


I 


(2sin»r  —  Dir 
(/'l— sin'r)*" 


s 


4«,^  R 

(A+p.>‘+(i,2+z.)*  • 


«  — 2r, 


(3) 


The  integral  in  expression  (3)  forms  the  sum  of  the  integrals, 

one  of  which  f— 7=^^=^  gives  us  V  =  during  the 

J  (K  1— *»s»ft*T)*  J  (K  1— **SJ||*T)» 

solution  by  parts. 


The  integral  V  is  solved  by  the  replacement  of  sUn  — t  ,  with  the 
subsequent  substitution  of  K(i  — -(i+U)*. 

As  a  result,  we  obtain: 

A  +  f  2-*» 

'*  ”  VT  12(1-  *•)  f^{p.  *)  -  f  1  - 

-  *.  *,  1  -  2Ti^)  T-X  ^  ^  ^ 

«fa2?  sin  2s 

Krrsnisrr/i’ 

where  E,  F  are  the  incomplete  elliptical  integrals  of  the  second  and 
first  types,  respectively,  taken  from  the  table  in  report  [6],  where 
Z'— (9)— f(v).  F(2.-H»±(p)— 2«/c±/(»);  K  is  the  total  integral.  The  situation  is 
analogous  for  E. 
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b)  With  Af  determined  from  (2)  with  the  replacement  of  cos 
by  — ,  the  component  of  induction  on  the  Z-axis  will  be: 


.  k*(0+l)-2  ,, 

2(1— A*)  A)  1  ~ 

_  D— 2  ^  /  iin  2? _ sin  :» 

-(!—**)  2  I  —  4»  sln‘ J  W-k-  sin'T 


c)  The  tangential  component  of  induction: 

!x,i(L'2^g,)k 


V /^-  -V-  ■■  _  ■  -J  , 

V  li^  I  —  A’  sia*  3  V  \  —  A*  sin*  » 


For  c^mductor  4,  expressions  (4)-(6)  are  valid  with  the  replacement 
of  (ZQ  +  L/2)  by  (zq  -  L/2)  and  th^  circuit  -i  (the  value' of  k  changes 
in  this  case) . 

The  point  p»=o 'is  a  singular  point  for  expressions  (4)-(6).  The 
following  expressions  of  the  induction  components  were  obtained  for 
this  point: 

^  .>*,<7?  (£.. 2  +  ?,)  sin  cos  .  ' 

2s(;?»  +  (i/2  +  r.)»]3'*  ’  I 

*''=«>  2s*r/?,4-(i,2  +  e,)*]J'2  ’  i 

g  .yff  (^/2  z,)  sin  sin  •,  j 

»(P.0)-  2RlR'+(Z./2  +  r,)>l^'*  ■  ' 


Components  of  Induction  from  the  Side  Part  of  the  Excitation  Coil 
(Circuit  3) .  Only 


-I/a  V  ^’•■*■‘’0  “  -  »•)  +  (*  -  - 


is  present  in  ( 1 ) . 


We  obtain  the  following  expressions: 


!»si7?sjn2*  _ ^^2  —  «« 

lA— ->+(^)' 


t  2  +  ?, 


1  — 'm*  sin*  J  4- 


[j _ ■ 

V A 


4 


H  <S<<>»4-ffcoi2P 

*•*  +  h)*  (i  —  ^ 


Li  — T, 


|/  l-m*>-«n*P  +  (  V+J,”  ' 


(9) 


/.,2-f  r. 


/ 


/  L  2  +  7, 

1  _„.5in‘p+  I 


where 


For  circuit  1,  formulae  (9)  are  valid  when  p  is  replaced  by  o 
and  with  the  current  -i.  The  components  of  the  t-th  excitation  coil 
are  determined  from  the  above  expressions  when  To  is  replaced  by 

i)+<Pol.  where  p  is  the  number  of  pairs  of  poles  in  the  excitation 
system  (with  the  angle  counted  in  the  positive  direction). 

Calculating  the  Main  Magnetic  Flux.  We  will  define  the  flow 
through  the  polar  surface  a2'bc4*da  (Fig,  1)  as  «r  -  i  Trf?  [7].  We  will 
express  ®  as  the  sum  of  the  fluxes  from  the  current  in  the  side  parts 
of  the  coil  as  ®<i  ,  and  in  the  front  parts  —  ®a  Based  on  the  sym¬ 
metry  of  the  magnetic  system,  the  sides  are  found  under  identical 
magnetic  conditions;  therefore,  da  and  be 

ip  m 

♦•  =  2  J  4,,,rf7,-2  ^  A,ydz,.  (10) 

-112  -il2 

where  A3^3»  is  the  vector  potential  from  the  current  in  circuit  3  on 
line  3'  of  the  armature  contour;  Ai^3*  -  from  the  current  in  circuit 
1  on  line  3 ' • 

The  vector  potential  from  the  current  in  circuit  m  (1  or  3)  on 
line  3' 

(„  -  L.'7)  +  -  L  g)» 

3'  -  -  4«  Lrt\  ^y  (f,  +  Lfiy  +  «»  ,.  (11) 

■ys 

where  z\  is  the  coordinate  of  the  current  point  on  line  3'  on  the 
Z-axis, 

3-  *  F*  +  p»’_  fo,  ^ 

^  3’  “  **  +  Po  -  COJ  (f ,  — 
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Solving  expression  (9),  we  obtain: 

***«  =  +  'W  (fs,  j'  —  li.i') 


where 


3'  -  2  -i- 4. 3-  -  2  -j/^ (~2^}  +  4. 3'  + 


l  —  L. 

,l-L  2 

+  -0—  In _ 


.  /  fi  —  Lv  '  +  i- ,  /  7r+Tv»  r 

1/  2  ]  /+!  ~2“+|/^  (~)  +4. 

[/m 


(12) 


/-i.  .  /  n-Lv- 

~^\/  [—)  +-.,3. 


i±i+ 


+  4.3' 


The  flux  from  the  front  parts 


*1 

2 


♦.  =■  2(1,  J  j,  —  ^4,  j,  )  <<f,. 


(13) 


The  value  of  the  vector  potential  is  found  by  solving  expression 
(2),  as  a  result  of  which  we  obtain  the  formula  for  calculating  the 
flow 


Where 


2  —  *5 

^  2*,  *.)  “  *,)  i  ~ 

_ ^ 

k,  *.)  ~  ^(1.  *,)1  ~  2*,  l^(»-  *.>  ~  -’ix.  *.)i  + 


*2-  fL-l\* 

(-T-  +<«  +  ?.)* 

\  “  y 


4?,/? 


45,^ 


+1?  +  ?.)' 


f,  — »/;>  is  the  angular  value  of  the  polar  spacing. 


(14) 


The  integral  of  expression  (14)  can  be  found  by  any  numerical 
method  (e.g.,  using  the  Simpson  rule  [8]). 


The  value  of  the  flow  from  the  t-th  pole  is  determined  from  the 
formulae  given  above,  with  •v*  replaced  by  (?„-.%>(:— the  sign  of 
expression  (14)  is  equal  to  (-1)^"^  in  this  case. 
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The  induction  component  Sj,  =  fU)  of  the  superconductive  excitation 
system  of  a  dipole  machine  was  calculated  using  the  above  formulae 
at  To-o  and  po»i4,:5- io-«  m  (for  the  surface  of  a  cryostat).  The 
results  of  the  calculation  are  compared  with  the  experimental  data  in 
Fig.  2. 

The  calculated  value  of  B,  at  point 
(>,=.3,9-io-*  m,  z  =  0  was  experimental¬ 
ly  verified.  The  results  agree,  within  the 
limits  of  the  precision  of  measurement  of 
induction  with  instrument  type  Ye-11. 

The  main  flux  of  the  machine  on  the 
calculated  diameter  of  the  armature  was 
determined  from  (12)  and  (111).  The  arc 
spacing  was  18°.  The  excitation  coil  was 
replaced  by  an  infinitely  thin  turn.  The 
divergence  of  the  experimental  value  of  the 
flux,  determined  from  the  machine's  EMF, 
from  the  calcualted  value  is  less  than  5%. 


CONCLUSIONS 

1.  The  proposed  formulae  can  be  used  in  the  electromagnetic 
calculation  of  an  electrical  machine  with  superconductive  excitation 
without  a  magnetic  circuit. 

2.  The  magnetic  flux  of  a  machine  with  an  armature  diameter  of 
100-150  mm  can  be  calculated  by  replacing  the  coil  cross  section  with 
one  infinitely  thin  conductor  in  the  geometric  center  of  the  cross 
section . 

3.  During  the  precise  calculation  of  induction  in  both  the  zone 
of  the  armature  winding,  and  in  the  points  close  to  the  excitation 
winding,  it  is  necessary  to  divide  the  coil  cross  section  into  several 


0  2,5  5,0  -iQ-lu 


Fig.  2.  Change  in  induc¬ 
tion  along  machine  at  po= 

=  14.75  cm,  <Pc-0- 
T  -  zone  of  fhont  part  of 
excitation  winding; 

_  -  experimental  curve; 

-  _  _  -  calculated  curve; 
.  -  computer  calcula¬ 
tion. 


elementary  cross  sections,  and  to  replace  each  by  an  infinitely  thin 
conductor.  The  breakdown  of  the  cross  section  into  more  than  16 
sections  virtually  does  not  change  the  result  of  the  calculation  at 
points  which  are  at  a  distance  equal  to  the  height  of  the  coil. 


4.  The  formulae  can  be  used  to  calculate  magnetic  systems  with 
saddle-shaped  MHD  converters. 
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